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At zero temperature the sublattice magnetization of tire quantum spin-1/2 Heisenberg antifer- 
romagnet on a body-centered cubic lattice with competing first and second neighbor exchange (Ji 
and J2) is investigated using the non- linear spin wave theory. The zero temperature phases of the 
model consist of a two sublattice Neel phase for small J2 (AFi) and a coUinear phase at large J2 
(AF2). We show that quartic corrections due to spin- wave interactions enhance the sublattice mag- 
netization in both the AFi and the AF2 phase. The magnetization corrections are prominent near 
the classical transition point of the model and in the J2 > Ji regime. The ground state energy with 
quartic interactions is also calculated. It is found that up to quartic corrections the first order phase 
transition (previously observed in this model) between the AFi and the AF2 phase survives. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee, 73.43.Nq 



I. INTRODUCTION 

In recent years thermodynamic properties of frus- 
trated quantum Heisenberg antiferromagnets have been 
of intense interest both theoretically and experimen- 
tally in condensed matter physicsi^ The phase diagram 
of the quantum spin-1/2 Heisenberg antiferromagnetic 
(AF) model on two-dimensional (2D) lattices with near- 
est neighbor (Ji) and next nearest neighbor interactions 

(J2) have been studied extensively by different meth- 
ods M^^^i^i^iiaai .12 ^a^^ie^is^^^ tljg sq^g^^g 

lattice with nearest neighbor (NN) exchange interaction 
only, the ground state is antiferromagnetically ordered 
at zero temperature. Addition of next nearest neighbor 
(NNN) interactions break the AF order. The competi- 
tion between the NN and NNN interactions for the square 
lattice is characterized by the frustration parameter p. It 
has been found that a quantum spin liquid phase exists 
between pic ~ 0.38 and p2c ~ 0.60. For p < pic the lat- 
tice is AF-ordered whereas for p > p2c a collinear phase 
emerges. In the collinear state the NN spins have a par- 
allel orientation in the vertical direction and antiparallel 
orientation in the horizontal direction or vice versa. 

Motivated by the results for the 2D lattices some work 
has been done by analytical and numerical techniques to 
understand the magnetic phase diagram of three dimen- 
sional (3D) latticesi^i^i^'2^i2^ Linear spin-wave theory, 
exact diagonalization, and linked-cluster series expan- 
sions (both at zero and finite temperature) have been 
utilized to study the 3D quantum spin-1/2 Heisenberg 
AF on a body-centered-cubic lattice (bcc) latticci^i^^ It 
has been found that the lattice does not have a quantum 
disordered phase and a first-order phase transition from 
the AF-phase (AFi) to lamellar state (AF2) occurs at 
Pc = 0.53 or J2/J1 ~ 0.705. The first-order nature of the 
phase transition from the AFi to the AF2 phase in the 
model is inferred from a kink in the ground state energy 
of the system. In one and 2D due to reduced phase space 



quantum fluctuations play an important role in deter- 
mining the quantum critical points of the system at low 
temperature. However, in 3D the phase space available 
is greater and quantum fluctuations play a lesser role. 
Hence, the absence of the quantum disordered phase for 
the BCC lattice. 

In this work, we study the 3D quantum spin-1/2 AF on 
a bcc lattice using the non-linear spin wave theory where 
we consider interactions between spin waves up to quar- 
tic terms in the Hamiltonian. We compute the effect of 
these higher order terms on the sublattice magnetization 
(see Fig[2]). The corrections to the magnetization become 
important as the classical transition point Pc=0.5 is ap- 
proached. Also, our calculations re-confirm the first order 
nature of the phase transition found in Refs. [2^l25l up to 
quartic interactions (see Fig [3]) . The paper is organized 
as follows. In Section|TT]we begin with a brief description 
of the properties of the bcc lattice relevant to our calcula- 
tions. We then set-up the Hamiltonian for the Heisenberg 
spin-1/2 AF on the bcc lattice. The classical ground state 
configurations of the model and the different phases are 
then discussed. Next we map the spin Hamiltonian to 
the Hamiltonian of interacting bosons and the non- linear 
spin-wave theory for the two phases are developed. The 
sublattice magnetizations and the ground state energies 
for the two phases are numerically calculated and the re- 
sults are plotted and discussed in Section IIIII Finally we 
summarize our results in Section [IVI 



II. FORMALISM 

Body-centered-cubic lattice consists of two interpen- 
etrating, identical simple cubic lattices, each of which 
consists of two interpenetrating, identical face-centered 
lattices. This makes the bcc lattice a 3D bi-bipartite 
cubic lattice. The basis vectors of the bcc lattice con- 
necting eight (zi = 8) nearest neighbors are (in units 
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(a) AFi phase 
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FIG. 1: AFi and AF2 ordered phases of the bcc lattice. In 
the AFi phase all A-sublattice spins point in the direction of 
an arbitrary unit vector while B-sublattice spins point in the 
opposite direction. For the AF2 phase there are two inter- 
penetrating Neel states each living on the initial sublattices 
A and B. 
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cos(fc^) cos(fcy) cos(fcz), 



(3) 



cos(2fcj;) + cos(2fcy) + cos(2fc2) /3. (4) 



where N is the number of sites on the lattice and p = 
Z2J2/Z1J1 is defined to be the parameter of frustration. 

At zero temperature, the classical ground state for the 
bcc lattice has two phases. In the limit of small p or 
J2 << Ji three isolated minima in energy, Eq/NJi — 
—45*^(1 — p) occur at the wave-vectors (tt, 0, 0), (0, tt, 0) 
and (0, 0,7r). They correspond to the classical two- 
sublattice Neel state (AFi phase) where all A-sublattice 
spins point in the direction of an arbitrary unit vector fi 
while _B-sublatticc spins point in the opposite direction 

n. 

In the other limit, for large p or J2 >> Ji, there is 
a single minimum in energy, Eq/NJi = —4:S^p at k = 
(7r/2, 7r/2, 7r/2). In this case the classical ground state 
consists of two interpenetrating Neel states (AF2 phase) 
each living on the initial sublattices A and B. The two 
phases are shown in Fig. [TJ 

The classical limit for the phase transition from AFi 
to AF2 for the 3D model on the bcc lattice is at the 
critical value Pc — Iz^jizx = 1/2 i.e. when J2/J1 — 2/3. 
This is similar to the spin- 1/2 J\ — J2 model on a 2D 
square lattice where the critical value of Pc — 1/2 or 
J2/J1 = 1/2. 



of simple cubic lattice spacing) ai = (1,1,-1), a2 = 
(1,-1,1), as = (—1,1,1) and the lattice vectors con- 
necting six (2,2 = 6) next-nearest neighbors are bi = 
(±2,0,0), ba = (0,±2,0) and bg = (0,0, ±2). On such 
a lattice the Hamiltonian for a spin-1/2 Heisenberg AF 
with first and second neighbor interactions is 

iJ=iji^S,.S, + ij2^S,.S,. (1) 

where Ji is the NN and J2 is the frustrating NNN ex- 
change constants. Both couplings are considered AF, i.e. 
Ji, J2 > 0. 

A. Classical ground state configurations 

The limit of infinite spin, S 00 corresponds to the 
classical Heisenberg model. We assume that the set of 
possible spin configurations of the system are described 
by Si = S'ue*'! '"', where u is a vector expressed in terms 
of an arbitrary orthonormal basis and q defines the rela- 
tive orientation of the spins on the lattice?^ The classical 
ground state energy of the system expressed as a function 
of the parameters Ji and J2 takes the form 

E^/NJi = i5'zi[7ik +P72k], (2) 



B. Non-linear spin wave theory 

The Hamiltonian in Eq.[T]can be mapped into an equiv- 
alent Hamiltonian of interacting bosons by transforming 
the spin operators to bosonic operators a, a' for A sub- 
lattice and b, for B sublattice using the well-known 
Holstein-Primakoff transformations^* 

Si. - V2S{1 ^) a., « V2Sal (l - ^) , 

^ V2^6](l-^), 5^^.«^/25(l-^)6„ 
5^^. = -S + b]b„ (5) 

In these transformations we have kept terms up to the 
order of 1/5. Next using the Fourier transforms 

k k 

the real space Hamiltonian is transformed to the k-space 
Hamiltonian. The reduced Brillouin zone contains N/2 k 
vectors as the unit cell is a magnetic supercell consisting 
of an ^-site and a i3-site. In the following two sections 
we study the cases J2 < Ji and J2 > Ji separately. 



3 



1. J2 < Ji: AFi phase 

In this phase the classical ground state is the two- 
sublattice Neel state (see Fig.[T]). For the NN interaction 
spins in A sublattice interacts with spins in B sublattice 
and vice- versa. On the other hand for the NNN exchange 
J2 connects spins on the same sublattice A with A and B 
with B. Substituting Eq. [5]in Eq. [1] expanding the rad- 
ical, and restricting to terms only up to the anharmonic 
quartic terms, we obtain the k-space Hamiltonian 



The quartic corrections to the ground state energy is 
calculated from the four-boson averages. In the leading 
order they are decoupled into the bilinear combinations 
fEqs. [T^- [H|) using Wick's theorem. The corresponding 
four boson terms are, 



{ala,b]bj) = it^ + t;^ {a\b]h.jhj) = 2uv, 
{a\aiaibj) = 2uw, {a\aio}jaj) = + in^ , 
{aiOjOjaj) — 2uw, (alaiaiOj) — 2uw. 



(18) 



i/ = i/(0)+i/(2)+i7(4). 



(6) 



The classical ground state energy H'^^^ and the quadratic 
terms i?^^-* are 



This yields the ground state energy correction from the 
quartic terms 



(7) 



(u + — p{u — w)' 



(19) 



k 



Summing all the corrections together the ground state 
energy takes the form 



+ Sok(ak^_k + a-k^k) , 

with the coefficients Aok and i?ok defined as 

^ok = l-p(l-72k), (9) 
Sok = 7ik. (10) 

The quartic terms in the Hamiltonian _ff are 
H^'^'' = - Ji ^ |^a|a,j6]6j -I- ^ (^aib^^bjbj + a\aiaibj + h.c. 



(8) E/NJ, = _izi5(5 + l)(l-p)-f izi^J^^L^k] 



—21 
2 ^ 



(u + v){l — u — v) — p{u — w)(l — u + w) 



(20) 



and the sublattice magnetization (Sa) at zero tempera- 
ture is given by 



{Sa 



2S\n ^ / 

k 



(21) 



M 
b] 



(11) 



Using Eqs. [151 HH ^^^d [T7] we numerically evaluate 
E/NJi and (Sa)- For the bcc lattice the k-sum is re- 
placed by an integral over the Brillouin zone^ 



These terms are evaluated by applying the Hartree-Fock 
decoupling process.— In the harmonic approximation the 
following Hartree-Fock averages are non-zero for the bcc- 
lattice Heisenberg AF: 



N 



E 



1 



"'O ^0 



(22) 



2iN ^ cjok 



{a,bj) 



2VN ^ 



Tlk^Ok 



Wok 



= {A^j)-{b\b,) = \[^Y. 

k 



72kA 



Ok 



i^Ok 



(12) 
(13) 
(14) 



where wok = \/^ok ^ ^ok- 

The contributions of the decoupled quartic terms to 
the harmonic Hamiltonian in Eq. [8] are to renormalize 
the values of ^ok and i?ok which are now 



2. J2 > Ji: AF2 phase 

The classical ground state for J2 > Ji corresponds to 
a four sublattice state where each of the A and B sub- 
lattice is itself antiferromagnetically ordered (see Fig.[T|). 
For the NN exchange there are four A — A, four B ~ B, 
and eight A — B type interactions between the sublat- 
tices. In case of NNN exchanges there are a total of 
twelve A — B type interactions. Adding all their contri- 
butions together up to the quadratic terms the harmonic 
Hamiltonian takes the same form as Eq. [8] with 



Ar. 



1 



u + V 



S 



p[l - 72k] 1 



S 



Sk = 7ik(l- 

Wk 



U + V 

S 



(15) 
(16) 
(17) 
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Aok 
-Bok 



(7ik + 2p), 
(71k + 2p72k) 



(23) 
(24) 
(25) 
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The quartic terms in the Hamiltonian for this case are 



H 



(4) 



— - Ji "^2, + ^ {^i^l^j^j + ajaittibj + h.c. 

+ -Ji ^ |^a-aia]aj - - (^aia^^ajaj + aJaiOiaj + h.c^ 
+ a ^ &j — J2 ^ |^ajai6]6j + - {aib^-bjh-j + ala^aibj 



h.c. 



(26) 



These terras are decoupled and evaluated in the same way 
as before. The renormalized values of the coefficients 
and i3k are 



Au = 



where 
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(27) 
(28) 

(29) 
(30) 



In Eqs. [27| and [55] u, v have the same form as in Eqs. [T^] 
~ [13] but they are evaluated with the coefficients Aok and 
Sok in Eqs. [Ml and [23 The quartic corrections to the 
ground state energy is 

^ iiVJiZi|(u- w)2 - {u + v)^ -p{u + vfy (31) 

Combining all these corrections the ground state energy 
is 



{v + w)(l — 2w — + w) 
+ 2p{u + v){l — u — v) . 



(32) 



The sublattice magnetization and ground state energy 
are then obtained numerically using Eqs. [nU [571 ESI 
and [31 



III. RESULTS 

In Fig. [1 we show the results for the sublattice mag- 
netization, (Sa), obtained numerically from Eg. 1211 for 
both AFi and AF2 phases with (dashed line) and with- 
out (solid line) quartic corrections. In the AFi ordered 
phase or the two sublattice Neel phase where A and B 



sublattice spins point in the opposite directions, sublat- 
tice magnetization decreases monotonically with increase 
in p till p Ri 0.5. The curve starts at « 0.44 for p = 
and ends at « 0.34 for p — 0.5. The gradual decrease in 
(Sa) is expected with increase in p as increasing strength 
of NNN interaction J2 aligns the spins antiferromagnet- 
ically along the horizontal and the vertical directions. 
The quartic corrections produce a change in the sub- 
lattice magnetization, (Sa), which becomes significant 
as one approaches the classical transition point Pc=0.5 
(see Fig. [5]). With quartic corrections the magnetization 
curve starts at « 0.44 for p = and ends at « 0.38 for 
p = 0.49. At p = (no frustration) there is no quartic 
corrections to {Sa}- This can be observed from Eqs. [TSJ- 
[T7][2T]as the correction factor {l — {u + v)/S) cancels out 
in Eq. [5TJ At the wave-vector k — (7r/2, 7r/2, 7r/2) spin- 
wave theory calculations become unstable (at Pc ~ 0.5) 
since the coefficient Ak becomes equal to _Bk- In the 
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FIG. 2; Sublattice magnetization, (5*^), is plotted versus p for 
AFi and AF2 ordered phases. In the AFi phase with increase 
in p the system aligns the spins antiferromagnetically along 
the horizontal and the vertical directions - thus decreasing the 
sublattice magnetization. In the AF2 phase {Sa) mostly stays 
the same and then shows a slight decrease (without quartic 
corrections) as p approaches the critical value Pc = 0.5 from 
above. However with the quartic corrections {Sa) remains 
almost constant at 0.43. In both cases quartic corrections 
to the Hamiltonian of the system enhance the magnetic order. 

AF2 ordered phase or the lamellar phase with two inter- 
penetrating Neel states, sublattice magnetization stays 
mostly flat except for a slight decrease (without quartic 
corrections) as p approaches the critical value pc from 
above. The curve starts at « 0.42 for p = 1 and ends at 
w 0.41 for p = 0.5. However with quartic corrections the 
curve has a very small upward turn. This upward curve 
has been observed in previous numerical works on this 
modeli^i^^ For the AF2 phase, quartic fluctuations pro- 
duce an overall enhancement of the magnetization over 
the high-p values (0.5 to 1). But for low-p (0 to 0.5), with 
increase in frustration quantum spin fluctuations play a 
dominant role as seen in Fig. [5] 
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FIG. 3: Ground state energy per site, E/NJi, is plotted as a 
function of the frustration parameter p = Z2J2/ z\J\ without 
(solid lines) and with (dashed lines) quartic corrections for 
both AFi (p < 0.5) and AF2 (p > 0.5) ordered phases. For 
the bcc-lattice zi = 8 and Z2 = 6. Spin wave theory becomes 
unstable at the classical transition point, i.e. p « 0.5. After 
extrapolation (not indicated in the figure above) we find that 
the two energies meet at p ~ 0.53 or J2/J1 ~ 0.705. The kink 
in the energy at this value of p indicates a first order quantum 
phase transition from AFi to AF2 phase. 

In Fig. [3] we plot the ground state energy per site, 
E/NJi, for the AFi and AF2 phases with and without 
quartic corrections as a function of the frustration param- 
eter p — Z2J1I z\3\. pc — 0.5 is the classical transition 
point where a phase transition from the AFi phase to the 
AF2 phase occur. The quadratic calculation agrees well 
with the results of [23l25l . The quartic corrections to the 
energy are shown by the dashed lines in Fig. [31 At p ~ Q 
calculated energy with the quartic correction is slightly 
lower than the energy calculated without the quartic in- 
teraction terms. This small decrease from the linear spin 
wave theory calculation is due to the ground state en- 
ergy correction which is negative (as seen in Eq. llOp from 
the quartic terms (self-energy Hartree diagrams). This 
trend for low-p continues till p ~ 0.38 after which the 
energy with quartic corrections become dominant. For 



large p, we find the energy with quartic corrections to 
be lower than the energy calculated without the quar- 
tic interactions in the interval « 0.70 — 1. In both the 
phases quantum spin fluctuations tend to maintain the 
magnetic order by lowering the ground state energies. 
As p approaches the critical value Pc from both phases, 
frustration increases causing the ground state energies to 
increase. Then 1/5' corrections due to spin fluctuations 
play a lesser role. As mentioned in the magnetization cal- 
culation our non-linear spin wave analysis becomes un- 
stable at the classical transition point Pc = 0.5. After 
extrapolation of the ground state energy curve from the 
AFi phase in the regime where non-linear spin wave the- 
ory breaks down we find that the energies from the two 
phases meet at p « 0.53.^^ The kink at this point signals 
a first-order phase transition occurs from AFi to AF2 
phase. 



IV. CONCLUSIONS 

In this work we have investigated the zero tempera- 
ture 1/S corrections to the sublattice magnetization and 
ground state energy of a spin-1/2 Heisenberg frustrated 
antiferromagnet on a bcc lattice using the framework of 
non-linear spin wave theory. We have found that 1/S cor- 
rections due to spin-wave interactions cause noticeable 
changes to the sublattice magnetization for both the two 
sublattice Neel phase (small NNN interaction J2) and the 
AF2 phase or the lamellar phase (large J2). As non-linear 
spin wave theory calculations become unstable close to 
the classical transition point we are unable to analyze 
the nature of phase transition using this method. We 
also confirm that up to quartic corrections the system 
undergoes a first-order phase transition as indicated by 
a kink in the energy calculation. 
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